Oscillating Fracture in Rubber 
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We have found an oscillating instability of fast-running cracks in thin rubber sheets. A well- 
defined transition from straight to oscillating cracks occurs as the amount of biaxial strain increases. 
Measurements of the amplitude and wavelength of the oscillation near the onset of this instability 
indicate that the instability is a Hopf bifurcation. 



PACS numbers: 62.20.Mk, 81.05. Lg, 83.60. Uv, 89.75.Kd 

When a balloon is pierced, a crack races around the 
circumference of the balloon, slicing the material into 
fragments. Surprisingly, instead of running straight the 
crack wiggles, leaving a wavy pattern on the fractured 
edges of the fragments. This effect has likely been ob- 
served many times, but we found only one documented 
case Q. We have constructed an experiment to study 
this phenomenon. 

Elastomers fill an important technological need, and 
understanding their properties is of commensurate im- 
portance. Aside from practical concerns, the oscillations 
of cracks in rubber may shed light on the broader issue 
of crack path selection. The field of dynamic fracture 
was spawned by theoretical attempts to understand the 
onset of crack branching in brittle materials B, an in- 
stability marked by the bifurcation of a straight-running 
crack. Despite the great advances in the understanding 
of dynamic fracture since then 0] , the crack path selec- 
tion problem remains unsolved. The oscillation of cracks 
in rubber may provide a new avenue to address this issue, 
as this is the first case we know where a rapidly moving 
crack spontaneously chooses a wavy path in a homoge- 
neous setting. 

There is a vast literature on the fracture of elastomers 
dating back a half century J4| . The applicability of the 
principles of fracture mechanics to elastomers has been 
established |5|~0], the time-dependent characteristics of 
elastomer fracture have been examined Mpl, an d crack 
speeds under varying degrees of biaxial strain have been 
investigated MiyJ. None of these studies, however, has 
explored the transition from a straight to a wavy crack 
path. Our aim here is to characterize this instability. 

Inflating a balloon introduces biaxial strain into the 
material and a pressure drop across the surface. The 
pressure difference plays no role so it is sufficient to ex- 
periment with flat sheets in biaxial tension. Our ap- 
paratus, inspired by a similar technique developed by 
Treloar |l2| , is shown Fig. n\ Most of our rubber samples 
were taken from a single roll of 0.18 mm thick natural 
latex sheet. Tabs were prepared on 32.5 cm x 12.7 cm 
sheets by cutting 1.2 cm long slots, perpendicular to the 



edge, 2.5 cm apart (Fig. 0(a)). The tips of these inci- 
sions were rounded by melting them with a soldering iron 
to prevent cracks from initiating at these points during 
loading of the sample. 
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FIG. 1. The experimental apparatus for straining rubber 
sheets along two axes, (a) A grid is drawn on the sample and 
clamps are attached to precut tabs along the sample's edges. 
The load is applied to the sample through the clamps, which 
are attached by wires to the rigid outer frame, (b) After the 
sheet has been slowly extended in the x- and y-directions, it is 
clamped by an inner rectangular frame. The sheet is distorted 
near its edges but not inside the inner frame. After clamping 
the sheet, it is punctured with a pin at the point marked x . 
Since the edges of the sheet are clamped, no energy flows into 
it during fracture. 



The experiment proceeds by gripping the tabs and in- 
crementally increasing the load simultaneously in both x- 
and y-directions, until the desired strain level is reached. 
The applied strain is on the order of 200%, uniform 
within 5%, and always chosen so that the strain in the y- 
direction, e y , is greater than the strain in the x-direction, 



e x . Strain is measured from the dilation of the grid; 
deviations from uniform strain are identified from the 
distortion of the grid and minimized by individually ad- 
justing each clamp. Once the desired strain level is at- 
tained, the rubber sheet is sandwiched between a pair of 
10 cm x 66 cm rectangular steel frames (Fig. 0(b)). The 
loading is then maintained entirely by the frames. 

Each run is initiated by pricking the sheet with a pin 
at the point marked x in Fig. Mb). The crack tip that 
forms is sharp and wedge-shaped, as shown in Fig. @. 
The crack travels down the centerline, the midpoint be- 
tween upper and lower edge of the frame. Depending on 
the initial strain conditions, the crack runs straight or 
oscillates about the centerline. The inner frame makes 
steady states possible because the energy stored per unit 
length inside the framed sample is constant ahead of the 
crack and also constant in its wake. So as the crack tip 
advances it consumes a fixed amount of energy per unit 
length of advancement. Indeed, we observe that after an 
initial transient period, an oscillating crack propagates 
with a wavelength and amplitude that are constant to 
within 10%. 



from any off-centerline starts, as shown in the inset of 
Fig. |(a). 

Rubber sheets with many different initial strain states 
were fractured, and the crack path was found to undergo 
a transition from straight to wavy with increasing biax- 
ial strain. Since the applied forces are purely tensile, 
the strain state is fully described by e x and e y . The re- 
sults of these runs are shown in the phase diagram in 
Fig. 0. The control parameter range was limited by ex- 
perimental difficulties found at the highest and the lowest 
values of e y . For e y > 2.6 it became impossible to com- 
plete a run because cracks would spontaneously form at 
the high stress point between the tabs. For e y < 1.4 it 
became impossible to distinguish between straight and 
wavy cracks because the wavelength became comparable 
to the length of the sample, as indicated by the trend 
in the inset of Fig. [|. Ambiguous points appear in 
Fig. ||(a) at high values of e y because different runs with 
the same initial conditions produced both straight and 
wavy cracks; ambiguous points at low values of e v result 
from the difficulty in discriminating between straight and 
wavy cracks. 




FIG. 2. Snapshot of the tip of a propagating crack. The 
crack is moving from left to right. The initial conditions were 
e y — 2.3 and e x = 1.8. 



Examples of the paths of a straight and a wavy crack 
are shown in Fig. Bl These curves were obtained by scan- 
ning the cracked sheet with a flatbed scanner and ap- 
plying an edge-finding algorithm to the resulting image. 
Transients typically dominate the first 15% of the crack 
length. During this transient regime a wavy crack's oscil- 
lations grow to saturation and a straight crack recovers 




= 2.0 (inset 

shows the initial kink that is sometimes observed) and (b) an 
oscillating crack with e x = 1.3, e y = 1.8 are shown in the 
final, unstrained state. In both cases, the crack initiated at 
the left edge and propagated to the right. Both samples were 
0.18 mm thick, held by a 10 cm x 66 cm frame. 

We measured the average wavelength (A) and the av- 
erage amplitude (A) of the wavy edge while holding 
€ y = 2.4 fixed and varying e x from 1.2 to 2.0. From the 
digitized curves we extracted the wavelength as the peak- 
to-peak distance in the x-direction and the amplitude as 
half the peak-to-valley distance in the y-direction, and 
averaged these quantities over multiple runs with identi- 
cal initial conditions. These data are plotted in Fig. 0. 

The amplitude grows as the square root of the control 
parameter with a critical value of x-strain, e x — 1.36. 
Furthermore, if we assume that the crack travels with a 



constant velocity v c in the x-direction, then the wave- 
length is equal to 2irv c /io, where to is the frequency at 
which the crack tip oscillates in the y-direction. Since 
at onset of the instability the wavelength is nonzero, the 
frequency at onset is nonzero. Hence, we conclude that 
the instability is a Hopf bifurcation. 
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FIG. 4. (a) A phase diagram showing the two observed 
states: oscillating cracks •, and straight cracks o, as a function 
of the initial strain state, e x and e y . The solid line shows our 
estimate of the phase boundary between straight and wavy 
cracks. Data were not taken in the region below the dashed 
line because there the principal direction of crack motion is 
across the width of the frame. (Frame size = 10 cm x 66 cm) 
(b) Wavelength versus e y , normalized by e x , for fixed e x = 1.4. 
The solid line is a fit of the data to A/(e y /e x — 1), which yields 
A = 0.58. 



The waveform of the crack path is nearly sinusoidal 
over a wide range of amplitudes (0.03 cm to 0.39 cm) 
and wavelengths (0.56 cm to 5.0 cm). Fig. @ illustrates 
this point with two waveforms taken from the runs used 
to construct Fig. u. These curves correspond to the first 
data point after the transition and the last data point 
in Fig. pi The waveform near the transition is almost a 
perfect sinusoid, consistent with the transition being a 
Hopf bifurcation. The waveform far from the transition 
shows sizeable deviations from a sinusoid and it is skewed 
in the direction of propagation. 

In addition to our quantitative results, we explored the 
possibility that the oscillation arises from out-of-plane 
vibrations, strain crystallization, and interaction of the 
crack tip with waves reflected from the boundary. In one 
experiment we reduced out-of-plane motion by sandwich- 
ing the rubber sheet snugly between two glass plates after 



the sheet was stretched to the loaded state. In another 
experiment we loaded the sheet as usual and then forced 
a cylindrical surface into the sheet so that the sheet was 
everywhere pressed into contact with the surface; thus 
the only out-of-plane motion possible is away from the 
surface. Neither experiment stopped the crack from os- 
cillating; hence out-of-plane motion is not the source of 
crack oscillations. 
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FIG. 5. Amplitude (A) • and wavelength (A) A are plot- 
ted across the transition from straight to oscillating cracks. 
e y = 2.4 (see Fig. M). The solid line through the amplitude 
data is a square-root fit. 




x/A. 
FIG. 6. Profiles of individual peaks from two different runs 
are scaled by their amplitude and wavelength and plotted 
together. Each run was performed under different strain con- 
ditions but for the same rubber thickness and frame size. A 
sine curve (solid line) is plotted for comparison. (•: e x — 1.4, 
€ V = 2.4, A = 0.03 cm, A = 0.75 cm; ♦: e x = 2.0, e y = 2.4, A 
= 0.13 cm, A = 1.02 cm). 



We also considered the effect of the known propensity 
of polymer chains in latex rubber to align when strained 
predominately in one direction p2[ . Measurements of 
the elastic response of our rubber samples indicate that 
the initial strains in our experiments were well below the 
threshold at which strain crystallization and finite chain 
length effects are significant 12 . Furthermore, by in- 
creasing the strain in both the x- and y-direction simul- 
taneously so that the strain is never much larger in one 
direction than the other, we avoided the conditions lead- 
ing to strain crystallization. We further considered the 
possibility that the high stress around the crack tip could 
locally crystallize the material as the tip passed through. 
However, the kinetics of the process is too slow. Assume 
that stresses are sufficiently large to crystallize the rub- 
ber out to some radius R from the crack tip. This region 
is traveling with a velocity v c and so the maximum time 
available for crystallization to occur is R/v c . The value 
of R must be less than the half- width of the frame, 5 cm, 
and the lowest crack velocity we measured was 40 m/s. 
This places an upper bound of 1.3 ms on time available 
for the rubber to crystallize. Comparing this result with 
the measured nucleation time of 60 ms for rubber at high 
strains [|13| , we find that there is insufficient time for the 
material around the crack tip to crystallize. A similar 
conclusion was reached by Lake and coworkers after find- 
ing similar crack growth rates for both a crystallizing and 
non-crystallizing natural rubber during high-speed frac- 
ture 0. 

In brittle materials, waves originating from the crack 
tip and reflecting from the boundaries are known to 
produce periodic markings such as those identified as 
Wallner lines K3]. If such a mechanism were active 
in rubber, one would expect an oscillatory wavelength 
A = (2hv c /c)/y/l — i>c/c 2 where h is the distance to the 
boundary. Given that the vertical boundary is 5 cm from 
the crack and that our measurements show v c /c < 0.5, 
it follows that A > 5.7 cm. Yet, the smallest wavelength 
observed was 0.85 cm (in the strained state). Thus waves 
reflecting from the boundary cannot account for the os- 
cillations of the crack. 

In conclusion, we have found a well-defined instability 
in the propagation direction of a crack in rubber sheets 
with biaxial strain. This is an unexpected result given the 
homogenous elastic field in which the crack propagates. 
We have shown that this instability can be characterized 
as a Hopf bifurcation. We have ruled out strain crystal- 
lization, out-of-plane motion and wave reflections from 
the boundary as possible mechanisms for the oscillation. 
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